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Abstract. The European Pulsar Timing Array (EPTA) is a multi-institutional, 
multi-telescope collaboration, with the goal of using high- precision pulsar timing 
to directly detect gravitational waves. In this article we discuss the EPTA member 
telescopes, current achieved timing precision, and near-future goals. We report 
a preliminary upper limit to the amplitude of a gravitational wave background. 
We also discuss the Large European Array for Pulsars, in which the five major 
European telescopes involved in pulsar timing will be combined to provide a 
coherent array that will give similar sensitivity to the Arecibo radio telescope, 
and larger sky coverage. 
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1. Introduction 



Pulsars are highly magneti zed neutron stars (NS), th e remnants of the supernovae of 
massive stellar progenitors dGoldl . 119681: IPaciml 1968). The majority of known pulsars 
rotate approximately once every second, and most are visible via their radio emission. 
It is generally accepted that radio pulsars radiate from above their magnetic poles; 
with each rotation of the NS, we obs erve a cross - section through the emission region 
that cuts our l i ne of sight at Earth ([Goldl . Il968t iPacinil . Il968t iGoldreich and Julian . 
19691 : ISturrockLll97lh . This gives rise to the "pulses" that lend pulsars their name. 
The intrinsic periodicity of most pu lsars is remarkably s t able, with typical period 



derivatives of ~ 10 ss (see e.g.. lLorimer and Krameii 120051 ) 



There is a population of so-called "recycled" pulsars with rotational periods 
of ~ 1 — 100 ms. They are the result of a mass- accretion episode from a stellar 
binary companion onto the NS. Th e accreted matter carries angular momentum, 
resulting in the "spin- up" of the NS ( Smarr and Blandfordl. 1976 : Alpar et al. . 1982 : 
Rad hakrishnan and Srinivasanl . Il982t ISrinivasan and van den Heuvell . Il982h . These 
millisecond pulsars (MSPs) are observed to be more stable as a population than 
the mainstream pulsars, with ro tational period derivatives of ~ 10~ 20 ss -1 (see 



e.g., Lorimer and Kramer! 120051) . They also show fewer instances of intrinsic, 



often low-frequency noise, corresponding to instabilities which affect the timing of 
mainstream pulsars. This is frequently referred to as "ti ming" or "red" n oise, and 



can be detrimental to long-term timing precis ion (e.g., lArzoumanian et all 11994 : 



Verbiest et al. . 2009t van Haasteren et all 2009). For these reasons, in addition to 



their 2 — 3 orders of magnitude shorter rotation rate, MSPs are thus the pulsars of 
choice when the highest timing precision is required, as is the case for the pursuit of 
gravitational wave (GW) detection. 

Through measurement of the pulse arrival times at Earth, we can continually 
improve our model ephemeris for the rotation of a pulsar. This is done by properly 
accounting for any deviations to the expected arrival times due to, for example, the 
motion of the Earth, the proper motion of the pulsar, or the binary parameters of the 
pulsar. This pulsar "timing" aims to successfully account for every rotation of the NS, 
which is referred to as "phase connection" ; this is what provides the high precision 
for which radio pulsar astronomy is known. The differences between the observed 
and expected times of arrival (TOAs) are referred to as "timing residuals" . An ideal 
timing model wi ll generally give residuals th at are normally distributed about a mean 
of zero (see e.g., Lorimer and Krameii 120051 ). 

Pulsar timing has been very successful for the measurement pulsar properties, and 
the study of the environments around NSs, for example. Perhaps its most celebrated 
results, however, are those that have tested general relativity (GR) and other theories 
of gravity, which we briefly outline in the next section. The potential that pulsar 
timing also holds for the direct detection of GWs is discussed in section 11.21 



1.1. Testing gravitational theories with pulsars 

The very first test of GR using pulsar observations involved timing of the first- 
discovered binary pulsar, PSR B1913+16. Through high-precision timing of this 
double-neutron-star system, the measured decay of its orbit was found to agree 
spectacularly well with the orbital e nergy loss due to quadrupolar GW radiation 
from that system, as predicted by GR ( Hulse and Taylor , 1975 : Taylor and Weisberd . 
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1989). More recently, observations of the double-pulsar system PSR J0737— 3039A/B, 
in which both NSs are observed as radio pu lsars, have provided the most stringent test 



to date of GR in the strong- field regime ([Kramer et all 120061 ) . In addition, several 



pulsars have been, and continue to be, used to test various aspects of GR and other 
theories of gravity. These include, for example, probing Strong Equivalence Principle 
violation by setting limits on phenomena predicted to exist by alternative theories of 
gravity. Among these are gravitational dipole radiation, secular evolution of binary 
parameters du e to preferred frame effects, and a variation in Newton ' s gravitational 



constant (e.g., Nordtyedtl. Il990t IPamour and Taylor . 1991 ; Wex . 2000t Arzoumanianl . 



2003; IVerbiest et all l2008t lLazaridis et all 12009). 



1.2. Direct detection of gravitational waves 

Although some of the tests described in the preceding subsection compare predictions 
of GR and other theories of gravity to pulsar observations in order to infer the 
existence of GW radiation as prescribed by GR, they fall short of being able to claim a 
direct detection of GWs. Pulsars can, however, be used to directly detect a stochastic 
GW background (GWB), pres umably produced by supermassive black hole binaries at 
the centers of distant galaxies ( Raiagopal and Roman! Il995t Jaffe and Backer . 200 



Damour and Vilenki 



2003|) 
snkinl . 



or pe rhaps by exotic phenomena such as cosmic strings (e.g 
2005). This is done by using an ensemble of pulsars as arms of a very large GW 
detector. The perturbing effect of a GWB passing by the Earth would be correlated 
amongst observations of all pulsars. This would manifest itself as a distinct signature 
in the timing residuals, once all model parameters particular to the pulsar system were 
properly accounted for. The cross-correlation of these residuals is predicted to show a 
distinct dependence on the separations between the various pulsars observed, and so 
a distribution of puls ars evenly s pread over as much of the sky a s possible is desired 
for such a study (e.g., Detweiler . 19791 Hellings and Downsi Il983l) . 

A distribution of sources in a "pulsar timing array" (PTA) thus provides a 
unique laboratory for performing this experiment. The PTA experiment is sensitive 
to GWs with frequencies of ~ 10 -9 Hz, since the limiting frequency corresponds to 
the timespan of the data sets used, typically 1 — 10 years in length. Such a detection 
will be complementary to th ose that will be made b y terrestri al and spac e -base d 
interferometers such as LIGO (|Abramovici et allll992h and LISA (|Danzmannl . l2000l ). 
which are most sensitive to wave frequencies of ~ 10 — 10 3 Hz and ~ 10~ 4 — 10~ 2 Hz, 
respectively. 

The expected correlated signal is very subtle, and its detection requires the 
highest-precision timing, with rms of timing residuals well below 1 /is. It is also 
crucial that the puls ars used for this undertaking are stable over long timescales (e.g., 



Verbiest et all l2009h . As stated in section [TJ most MSPs fit this description. Ideally, 
observations of 20 — 40 MSPs are required to be observed over 5 — 10 years, with 
consistent ~ 100 ns rms for the resulting timing residuals, in order to be able to 
significantly distin guish a GWB signal in the measured cross-correlation spectrum 



(|jenet et all [2005). 



The direct detection of GWs using a PTA is certainly a very challenging task. 
There now exist three major research groups dedicated to this effort: The Parkes 
Pulsar Timing Array (PPTA) in Australia, which uses the Parkes radio telescope 
( Verbiest et all l201Clh ; the North American Nanohertz Observatory for Gravitational 



Waves (NANOGrav) collaboration in North America, which takes its data from the 
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Table J 


. Comparison of telescopes taking part in PTA projects. 


Telescope 


Diameter 


Aperture 


System 


Allocated 


Declination 




M 


efficiency a 


temperature 21 (K) 


time (h/mo) 


range (°) 


EPTA member telescopes 


Effclsberg 


100 


0.54 


24 


24 


> -30 


Lovell 


76.2 


0.55 


30 


48 


> -35 


Nangay 


94 b 


0.48 


35 


250 


> -39 


Sardinia 


64 


0.60 


25 


30 


> -46 


WSRT 


96 b 


0.54 


29 


32 


> -30 


LEAP C 


194 


0.54 


30 


24 


> -39 d 



Other PTA telescopes 



Arecibo 


305 


0.50 


30 


8 


[-1.+38] 


GBT 


100 


0.70 


20 


18 


> -46 


Parkes 


64 


0.60 


25 


100 


< +26 



a Values are for 1.4 GHz. 
b Circular-dish equivalent diameter. 

c First light expected in 2010; shown are projected values. 

d Thc limiting declination range from combining at least two EPTA member telescopes. 



Arecibo and Green Bank telescopes ( Jenet et al. . 20091) ; and the European Pulsar 



Timing Array (EPTA). These groups have begun to develop a partnership, with the 
goal of sharing resources in an International Pulsar Timing Ar ray (IPTA) effort, w hich 



will be necessary to achieve the objective of GW detection ( Hobbs et al. . 2010f) . In 



the following section, we focus on the instruments used by the EPTA and its current 
progress; in section[3l we describe the Large European Array for Pulsars (LEAP), the 
next step for the EPTA toward GW detection; finally, we conclude in section |4] with 
a short summary. 



2. The European Pulsar Timing Array 

The EPTA is a multi-national European collaboration of pulsar astronomers with the 
common aim of performing and improving high-precision pulsar timing, principally 
for the purpose of detecting GWs. It is also the goal of the EPTA to discover 
new candidate PTA pulsars, and as such, EPTA members are involved in various 
international pulsar search programs. The EPTA consists of members from various 
institutions in France, Germany, Italy, the Netherlands, and the United Kingdom, 
with continued ties to research groups in other countries around the world. 



2.1. Telescopes 

The EPTA is unique among PTA collaborations, in that it employs the data from five 
large-diameter radio telescopes, all of which have ongoing pulsar observing programs, 
including those focused specifically on collecting MSP data for the purposes of a 
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PTA. In what follows, we describe the key features of each of these observatories. A 
summary of properties of the five telescopes described above, and a comparison with 
the other PTA telescopes, can be found in Tab le [J For a d i rect c omparison of the 



timing capabilities of the three PTA groups, see lHobbs et al. I (|2010h . 



Effelsberg. The Effelsberg radio telescope, administered by the Max-Planck- 
Institut fur Radioastronomie, is located in Effelsberg, Germany. At 100 meters 
in diameter, it is the world's second-largest fully steerable radio telescope, and 
can perform observations at frequencies between 0.365 and 95.5 GHz. For pulsar 
observations, it hosts a digital filterbank (DFB) system that can take data over 1 
GHz of observing bandwidth. It is also capable of performing coherent dedispersion 
over 112 MHz of bandwidth, and provides full polarization data, with plans to expand 
this observing window in the near future. 

The Effelsberg telescope has recently installed a new seven-beam focal plane L- 
band (~ 1.4 GHz) receiver. This will most notably be used in an upcoming pulsar 
search campaign, scheduled to begin in December 2009, in which about 500 "normal" 
pulsars, and 100 MSPs are expected to be found. 



Lovell. The 76-meter Lovell radio telescope is located at the Jodrell Bank 
Observatory in the United Kingdom. Its operation is overseen by the Jodrell Bank 
Centre for Astrophysics at the University of Manchester. Observing primarily at 0.35 
and 1.4 GHz, it is a key member of the MERLIN and European VLBI networks. 

For pulsar research, it has a DFB backend similar to that at the Effelsberg 
telescope, and is also capable of wide-bandwidth observations. The addition of a 
coherent dedispersion instrument is also being implemented, in order to enable even 
higher-precision data to be taken for pulsar timing observations. 



Nangay. Located in central France, the 94-meter-diameter-equivalent Nangay radio 
telescope is the fourth-largest single-dish telescope in the world, and observes at 
frequencies ranging between 1.1 and 3.5 GHz. It currently houses the Berkeley- 
Orleans-Nangay pulsar backend, which is currently designed to coherently dedisperse 
pulsar data over 128 MHz of bandwidth; this will be expanded to cover 400 MHz 
of observing bandwidth in the near future (e.g., ICognard et all [20091) . There is also 
a DFB mode that can perform observations over a similar bandwidth, and is used 
primarily for search observations. 

Although its meridian-style design restricts the telescope to track a given source 
for about 1 hour, pulsar observations are heavily scheduled at the Nangay telescope, 
amounting to an average of 250 hours of telescope time per month. This, along with 
its state-of-the-art instrumentation, makes this telescope very well-equipped for long- 
term pulsar timing observations. 



Westerbork Synthesis Radio Telescope (WSRT). Operated by ASTRON in 
the Netherlands, The WSRT is a 14 x 25-meter telescope array, whose signal is 
combined in phase to form the equivalent of a single 96-meter radio telescope. The 
WSRT can currently observe at frequencies between 0.12 and 8.7 GHz, of which 
pulsar timing data is regularly collected at 0.35, 1.4, a nd 2.3 GHz. The WSRT i s 
equipped with the PuMa-II pulsar backend instrument (|Karuppusamv et al. L l2008h . 
which performs coherent dedispersion of the incoming pulsar signal over a bandwidth 
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Table 2. Sample of MSPs currently monitored by the EPTA. 




Pulsar 


Total data 


rms of 




name 


time span (yr) 


residuals (/xs) 




PSR B1640+22 


12 


1.6 




PSR J1713+0747 


11 


0.73 




PSR J1744-1134 


10 


0.55 




PSR B1855+09 


23 


1.7 




PSR J1909-3744 


4 


0.11 




PSR J1918-0642 


7 


2.2 





of up to 160 MHz. A major upgrade is in the works for outfitting each of the 
individual tele scopes with a focal-pl ane array receiver, called APERTIF (APERture 



Tile In Focus; IVerheijen et all 120081 ). This will effectively give the WSRT a factor of 
25 increase in field of view. It will reduce the observable frequency window to range 
between 1.0 — 1.7 GHz, but will allow for 300 MHz of continuous observing bandwidth. 

Sardinia. In Italy, the National Institute of Astrophysics (INAF; Istituto Nazionale 
di Astrofisica), is constructing the Sardi nia Radio Telescope (SRT), a fully-steerablc 
dish that will be 64 meters in diameter ( Tofani et all 120081 ) . Most useful for pulsar 



research will be its planned dual-receiver system that will be capable of performing 
simultaneous observations at P-band (305 — 425 MHz) and L-band (1.3 — 1.8 GHz). 
At the present time it is planned that the SRT will house a DFB backend for pulsar 
observations. First light for the SRT is planned for 2010. 

The EPTA is distinct among PTA collaborations, in that it continues to 
provide valuable experience in combining data from several observatories. This will 
demonstr ate how to effective ly and simultaneously deal with these varied data sets 
(see, e.g., Janssen et all [20081) . This will be of great importance for a future combined 



effort of an IPTA collaboration that will combine the data and resources of the 
NANOGrav and PPTA groups with those of the EPTA. 

2.2. Current progress 

All the currently operational telescopes of the EPTA have observing programs that 
involve the regular monitoring of MSPs. The quality of the data, and timing precision 
obtained, is comparable with those of the other PTA collaborations. Table [2] shows a 
sample of six MSPs currently observed by the EPTA, their respective data set time 
spans, and timing residual rms values. 

There is also an active effort to model the observed timing residuals 
into a measurement of a GWB amplitude that is constrained by the EPTA 
data. This involves a Bayesian statistical analysis, employing a Markov Chain 
Monte Carlo algorithm to sample the likelihood function of the GWB amplitude 
( van Haasteren et al. . 2009). At the present time, the timing precision of the EPTA 



pulsars is insufficient to achieve a significant GWB detection. However, using this 
method, we obtain a preliminary 95% upper limit to the GWB amplitude of 1.9 x 10~ 14 
y r -i/2. likelihood function derived from this analysis is shown in Figure Q] For 
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1 2 3 

GWB amplitude (10" 



4 

1/2) 



Figure 1. Likelihood func- 
tion of the GWB amplitude 
derived from EPTA pulsar 
timing data, normalized so 
that the peak of the function 
has a value of 1.0. The dashed 
vertical line shows the 95% 
upper limit to the GW am- 



plitude of 1.9 X 10" 



-1/2 



corresponding to a GW en- 
ergy density /iq^qw = 2.1 x 
10 — 7 for a wave frequency of 
1 yr" 1 . 



a wave with a frequency of of 1 yr -1 , this corresponds to a GW energy density 
/i§J1gw = 2.1 x 10~ 7 (van Haasteren et al. 2010, in prep.). With more data, this 
limit will continue to become further constrained. 



3. LEAP: The Large European Array for Pulsars 

The existing collaboration of EPTA member institutions, together with the close 
proximity of 100-meter-class radio telescopes from each other, as well as the already 
well-established pulsar observing programs at these telescopes, presents a unique 
opportunity to combine the signals from the major European observatories in a 
coherent array. Such an endeavor will provide a significant boost to the quality of 
data used by the EPTA, enhancing its potential for detecting a GWB. This is the 
primary objective of the LEAP project, which is expected to see first light by the 
end of 2010. This section provides a short introduction to this project; for a more 
in-depth reference on LEAP, its capabilities, anticipated results, and contribution to 
PTA studies, see Stappers et al (2010, in prep.). 

LEAP will utilize techniques adapted from exis ting long-baseline arrays such as 
VLBI (e.g., iDeller et all . l2007t iBrisken et all , |2010|) and MERLIN. However, unlike 



the imaging data that are obtained from those collaborations, LEAP will be designed 
specifically for taking pulsar timing data. The goal is to produce a time series from 
all the telescopes added coherently in phase, significantly increasing the sensitivity 
of pulsar observations. This will require techniques beyond the current state of 
the art in data-handling and processing. The "tied array" formed by the EPTA 
member telescopes will correspond to the equivalent of a single dish of approximately 
194 meters in diameter or 29 500 m 2 of collecting area, comparable to that of the 
illuminated Arecibo radio telescope. This will substantially improve the timing results 
of the current timing array pulsars. Perhaps more importantly, the sensitivity of LEAP 
will allow for the precise study of MSPs that have otherwise been too weak for PTA 
study with other telescopes, including the individual member observatories of the 
EPTA. An important advantage that LEAP will have over the Arecibo telescope will 
be its much larger sky coverage. This will allow access to a greater number of currently- 
known pulsars — as well as those that will be discovered in the future — that can be 
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Figure 2. Comparison of the sky coverage limits of the Arecibo (dash-dot-dot-dot 
lines) and Parkes (strong dashed lines) telescopes with the future LEAP telescope 
array (strong solid line). Plotted are known MSPs, with those visible from Europe 
(and thus LEAP) labelled as stars. The central vertical and horizontal solid lines 
represent 0° in 1 and b, respectively. Divisions of 60° in 1 and 30° in b arc 
represented by weak dashed lines. 



used as part of a PTA (see Table [T] for a comparison between telescopes) . Figure [2] 
shows the distribution of known MSPs on the sky, and compares the ranges of sky 
visible from Arecibo, Parkes, and Europe (and thus LEAP). This will make LEAP the 
largest fully-steerable telescope in the world, and among the most sensitive, producing 
ideal data for high-precision puls ar timing and a major step forward in reaching PTA 



goals (e.g., Stappers et al. . 2009). 



The initial LEAP observing strategy will be to observe 20 — 30 MSPs over a 24- 
hour period on an approximately monthly basis for 3 years. This will be done over 100 
MHz of bandwidth, at an observing frequency of 1.4 GHz. The Nyquist-sampled raw 
voltage data will be recorded to hard disk, both from the MSPs themselves, and from 
calibrator sources. The latter of these data sets will provide the phase information that 
is needed for the multi-telescope tied-array to be achieved. The final data-combining 
step will be performed at the University of Manchester, where a coherent sum of 
the data in the directions of the observed pulsars will be calculated using their on- 
site supercomputer. The resulting high resolution, high sensitivity time-series data 
will then be processed by more familiar software to generate coherently-dedispersed, 
integrated pulse profiles. The final step will be to calculate the pulse TOAs, which 
are ultimately used for the high-precision timing program to detect GWs. 

However, this project will have its share of challenges; principal among these will 
be achieving a coherent sum of the raw recorded time-series data from all participating 
telescopes. Similar procedures are used for doing this at the WSRT; however, LEAP 
will be composed of dishes at distances over two orders of magnitude further away 
from each other than those of the WSRT, each having very different properties and 
observing systems. Furthermore, the tied-array will be formed in software rather than 
hardware. While this will provide greater flexibility than a system that combines 
the raw data in hardware, software-combining techniques are not yet sufficiently well- 
established for this type of data set. There will also be a challenge in data management 
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that must be overcome if LEAP is to be successful, as the volume of data that will be 
produced in a single observing session is expected to be in the hundreds of terabytes. 
Other issues with which to contend will be the difficult, yet crucial tasks of ionospheric 
calibration, characterization of radio frequency interference, and identification of the 
various other systematic effects on the data that will likely arise. 

All these potential difficulties will undoubtedly be encountered by future radio 
astronomy projects such as the Square Kilometer Array (SKA). Thus, in addition 
to advancing the PTA effort to directly detect GW radiation, LEAP will provide a 
much-needed test-bed for the future of radio observing techniques and precise pulsar 
timing with combined-signal telescope arrays. 

4. Summary 

We have outlined in this article some of the challenges that are faced by the pulsar 
astronomy community in detecting a stochastic GWB, and have described the role of 
the EPTA within this international effort. The EPTA, with its access to five large 
radio telescopes, already has the ability to provide high-quality pulsar timing data that 
is comparable to the other existing PTA groups. It is also in the process of moving a 
step further with the construction of the LEAP telescope array, which will coherently 
combine the time-series data from the five major European observatories. This will 
form a telescope that will be equivalent in sensitivity to the Arecibo radio telescope, 
but with much larger sky coverage. This project will not only provide very high- 
precision pulsar data, but will also be a vital testing ground for future multi-telescope 
arrays, most important of which is the SKA. 

Acknowledgements 

We are very grateful to all staff at the Effelsberg, Jodrell Bank, Nangay, and 
Westerbork radio telescopes for their help with observations used in this work. 
J. W. T. Hessels is an Netherlands Foundation for Scientific Research (NWO) Veni 
Fellow. K. Lazaridis was supported for this research through a stipend from the 
International Max Planck Research School (IMPRS) for Astronomy and Astrophysics 
at the Universities of Bonn and Cologne. Y. Levin and R. van Haasteren are supported 
by the NWO through Vidi Grant 639.042.607. M. Purvcr is supported by a grant 
from the Science and Technology Facilities Council (STFC). S. Sanidas is funded by 
a STFC Doctoral Training Account (DTA) scholarship. Part of this work is based on 
observations with the 100-m telescope of the Max-Planck-Institut fur Radioastronomic 
(MPIfR) at Effelsberg. Access to the Lovell telescope is supported through an STFC 
rolling grant. The Nangay radio telescope is part of the Paris Observatory, associated 
with the Centre National dc la Recherche Scientifique (CNRS), and partially supported 
by the Region Centre in France. The Sardinia Radio Telescope is a project of the 
Italian National Institute for Astrophysics (INAF). The Westerbork Synthesis Radio 
Telescope is operated by the Netherlands Foundation for Research in Astronomy 
(ASTRON) with support from the NWO. The LEAP project is funded through an 
European Reseach Council (ERC) grant to M. Kramer. 

References 

Abramovici A, et al 1992 Science 256, 325-333. 



REFERENCES 



10 



Alpar M A, Cheng A F, Ruderman M A and Shaham J 1982 Nature 300, 728-730. 

Arzoumanian Z 2003 in M Bailes, D. J Nice and S Thorsett, eds, 'Radio Pulsars', Astronomical 
Society of the Pacific, San Francisco, Vol. 302, pp. 69—74. 

Arzoumanian Z, Nice D J, Taylor J H and Thorsett S E 1994 ApJ 422, 671-680. 
Brisken W F, Macquart J, Gao J J, Rickett B J, Coles W A, Deller A T, Tingay S J and West C J 
2010 ApJ 708, 232-243. 

Cognard I, Theureau G, Desvignes G and Ferdman R 2009 in 'Windows on the Universe', XXI 
Rencontres de Blois, France, 21-26 June 2009 l|arXiv:0911.1612|l 

Damour T and Taylor J H 1991 ApJ 366, 501-511. 

Damour T and Vilenkin A 2005 Phys. Rev. D 71(6), 063510. 

Danzmann K 2000 Adv. Space Res. 25, 1129-1136. 

Deller A T, Tingay S J, Bailes M and West C 2007 PASP 119, 318-336. 

Detweiler S 1979 ApJ 234, 1100. 

Gold T 1968 Nature 218, 731-732. 

Goldreich P and Julian W H 1969 ApJ 157, 869-880. 

Hellings R W and Downs G S 1983 ApJ 265, L39. 

Hobbs G, et al. 2010 Class. Quantum Grav. 27 ( arXiv:0911.5206|l 

Hulse R A and Taylor J H 1975 ApJ 195, L51-L53. 

Jaffe A H and Backer D C 2003 ApJ 583, 616-631. 

Janssen G H, Stappers B W, Kramer M, Nice D J, Jessner A, Cognard I and Purver M B 2008 A&A 
490, 753-761. 

Jenet F A, Hobbs G B, Lee K J and Manchester R N 2005 ApJ 625, L123-L126. 
Jenet F, et al. 2009. larXiv:0909. 10581 

Karuppusamy R, Stappers B and van Straten W 2008 PASP 120, 191-202. 
Kramer M, et al 2006 Science 314, 97-102. 
Lazaridis K, et al 2009 MNRAS 400, 805-814. 

Lorimer D R and Kramer M 2005 Handbook of Pulsar Astronomy Cambridge University Press. 
Nordtvedt K 1990 Phys. Rev. Lett. 65, 953-956. 
Pacini F 1968 Nature 219, 145-146. 

Radhakrishnan V and Srinivasan G 1982 Curr. Sci. 51, 1096. 
Rajagopal M and Romani R W 1995 ApJ 446, 543-549. 
Smarr L L and Blandford R 1976 ApJ 207, 574-588. 
Srinivasan G and van den Heuvel E P J 1982 A&A 108, 143-147. 

Stappers B, Vlemmings W and Kramer M 2009 in 'Proceedings of the 8th International e-VLBI 
Workshop', 22-26 June 2009. Madrid, Spain p. 20. 

Sturrock P A 1971 ApJ 164, 529-556. 

Taylor J H and Weisberg J M 1989 ApJ 345, 434-450. 

Tofani G, et al 2008 in 'Society of Photo-Optical Instrumentation Engineers (SPIE) Conference 
Series' Vol. 7012, p. 70120F-12. 

van Haasteren R, Levin Y, McDonald P and Lu T 2009 MNRAS 395, 1005-1014. 

Verbiest J P W, et al. 2008 ApJ 679, 675-680. 

Vcrbicst J P W, et al. 2009 MNRAS 400, 951-968. 

Verbiest J P W, et al. 2010 Class. Quantum Grav. 27 HarXiv:0912 .2692 ) 

Verheijen M A W, Oosterloo T A, van Cappellen W A, Bakker L, Ivashina M V and van der Hulst 
J M 2008 in R. Minchin & E. Momjian, ed., 'The Evolution of Galaxies Through the Neutral 
Hydrogen Window' Vol. 1035 of American Institute of Physics Conference Series pp. 265-271. 

Wex N 2000 in M Kramer, N Wex and R Wielebinski, eds, 'Pulsar Astronomy - 2000 and Beyond, 
IAU Colloquium 177' Astronomical Society of the Pacific San Francisco, Vol. 202, pp. 113-116. 



